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Abstract 

In previous work, we introduced a new framework to treat large scale backreaction effects due to 
small scale inhomogeneities in general relativity. We considered one-parameter families of space- 
times for which such backreaction effects can occur, and we proved that, provided the weak energy 
condition on matter is satisfied, the leading effect of small scale inhomogeneities on large scale 
dynamics is to produce a traceless effective stress-energy tensor that itself satisfies the weak energy 
condition. In this work, we illustrate the nature of our framework by providing two explicit exam- 
ples of one-parameter families with backreaction. The first, based on previous work of Berger, is a 
family of polarized vacuum Gowdy spacetimes on a torus, which satisfies all of the assumptions of 
our framework. As the parameter approaches its limiting value, the metric uniformly approaches a 
smooth background metric, but spacetime derivatives of the deviation of the metric from the back- 
ground metric do not converge uniformly to zero. The limiting metric has nontrivial backreaction 
from the small scale inhomogeneities, with an effective stress-energy that is traceless and satisfies 
the weak energy condition, in accord with our theorems. Our second one-parameter family consists 
of metrics which have a uniform Friedmann-Lemaitre-Robertson- Walker limit. This family satisfies 
all of our assumptions with the exception of the weak energy condition for matter. In this case, 
the limiting metric has an effective stress-energy tensor which is not traceless. We emphasize the 
importance of imposing energy conditions on matter in studies of backreaction. 
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I. INTRODUCTION 

The standard model of cosmology is based on a hierarchy of length scales. At the largest 
scales in cosmology, it is generally believed that the universe is well-described by a met- 
ric with Friedmann-Lemaitre- Robertson- Walker (FLRW) symmetry that satisfies Einstein's 
equation with source given by the averaged stress-energy of matter. At scales somewhat 
smaller than this, the deviations from the FLRW solution are believed to be well described 
by linear perturbation theory. Finally, below a certain scale, nonlinear effects become im- 
portant, and numerical iV-body simulations are normally used to make predictions. The 
above beliefs underlie the standard cosmological model, which has proven to be extremely 
successful to date. 

Despite the success of the standard cosmological model, objections have been raised that 
nonlinear effects due to the small scale structure are not being fully taken into account 
(see, e.g., (l|). The Friedmann equations neglect averages of terms nonlinear in small scale 
perturbations, which are present in the Einstein equation. Since nonlinearities are important 
for small scale dynamics, how do we know that their averaged effects cannot contribute at 
large scales? Might it even be possible that such nonlinearities could produce effects on 
large scale dynamics that mimic the effect of a cosmological constant? 

At small scales, where nonlinear dynamics are important, fractional density perturbations 
can be much larger than 1. This formed the basis of Ellis's original argument |l| that 
neglecting nonlinear terms in the Einstein equation could be problematic. However, Ishibashi 
and Wald 2] argued that, even though perturbations in the matter stress-energy can be large, 
perturbations in the metric (as opposed to derivatives of the metric) should nevertheless be 
small, provided that appropriate gauge choices are made. They argued that, except in the 
immediate vicinity of strong field objects such as black holes and neutron stars, our universe 
should be well described on all scales by a Newtonianly perturbed FLRW metric, and that 
the terms in Einstein's equation that are nonlinear in the deviation of the metric from an 
FLRW metric should be negligible. 

However, it is clear that, in general, non-Newtonian circumstances, even when the devia- 
tion of the metric from a metric with FLRW symmetry is small, the nonlinear terms in the 
deviation of the metric from FLRW symmetry can affect the large scale dynamics. Indeed, if 
the universe were filled with low amplitude, high frequency gravitational radiation, then the 



metric could be uniformly close to a metric with FLRW symmetry, but this FLRW metric 
would behave dynamically as though the universe were filled with a P = |p fluid. This 
dynamical behavior occurs as a consequence of the nonlinear terms in Einstein's equation. 
Thus, an interesting issue in general relativity is to quantify the degree to which small scale 
spacetime inhomogeneities can, through nonlinear interactions, affect large scale dynamics. 
Prior to our work [3[, a number of averaging formalisms had been developed to try to 
address this issue. The most widely studied approach was initiated by Buchert |4|, [5]. Here, 
one works in comoving, synchronous coordinates, so that the metric takes the form 

ds 2 = —dt 2 + q i j{t 1 x)dx l 'dx 3 '. (1-1) 

One can then define averages of scalars over (compact, comoving regions of) constant-t 
hypersurfaces. In particular, Buchert defines an averaged "scale factor" as the cube root 
of the volume of one of these regions. Equations of motion for the averaged quantities are 
obtained by averaging (components of) Einstein's equation. In particular, by averaging the 
Hamiltonian constraint and the Raychaudhuri equation, analogs of the Friedmann equations 
can be derived in this framework. 

There are two independent serious difficulties with this framework. First, it is far from 
clear how to physically interpret the averaged quantities defined in this framework. For 
example, if one were to set up synchronous coordinates of the form (II .ip in our solar system 
based upon a family of geodesies that are initially (nearly) "at rest" with respect to the 
sun, then the geodesies near the planets would immediately begin shearing; the nonlinear 
effects of shear would soon result in significant convergence and eventually there would be 
caustics. The dynamics of the "scale factor" (i.e., the volume of a co- moving region) would 
be quite complicated and would show large effects produced by the planets. One might then 
be tempted to conclude that nonlinear effects produced by the planets result in large effects 
on the spacetime metric of the solar system. Obviously, this is not the case: The "large 
effects" are gauge artifacts resulting from the use of synchronous coordinates rather than 
coordinates in which the metric is (nearly) stationary, and quantities like the "scale factor" 
do not have any physical interpretation in terms of the observations of the metric made by 
(nearly) stationary observers. 

A second serious difficulty of this approach is that the "Friedmann equations" derived in 
this framework are not closed. They contain a scalar quantity Qt>, arising from averaging 



nonlinear terms in the Einstein equation, which is called the kinematical backreaction scalar. 
The evolution of Qv is undetermined in the framework. This is because the framework 
considers only the spatial average of two components of Einstein's equation, and thus is 
discarding much of the content of Einstein's equation. 

Another approach to backreaction, known as macroscopic gravity, was initiated by Za- 
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laletdinov [6|. Rather than restricting to averaging scalar quantities on spatial slices, in 
macroscopic gravity, an averaging operator (•) is introduced which allows one to define local 
spacetime averages of tensor quantities. The idea is to derive an equation for the average 
{gab) of the actual spacetime metric g a t,. Backreaction in this framework occurs through the 
"connection correlation tensor" , 

z\ c e f — (r a &[ c r | e |j]) — (r a b j c )(r | e |^), (1.2) 

where r a fec is the Christoffel connection. Various properties are derived for Z a b c def, but, as 
with <2x>, the full content of Einstein's equation is not used. 

Because of the presence of undetermined quantities in the averaged frameworks of Buchert 
and Zalaletdinov, a common approach to "solve" the equations has been to assume a par- 
ticular form for the backreaction tensors, and from that starting point, derive results for 
the large scale behavior 7|, |8| . However, no argument is given that the assumed form of the 
backreaction tensors can actually arise from averaging a physically reasonable small scale 
matter distribution. For example, in [8|, Buchert and Obadia assume that Qt> is equivalent 
to a scalar field (called the morphon field [9]), with a particular potential. With a suitable 
choice of potential, they find that inflation can occur in vacuum spacetimes due to the pres- 
ence of inhomogeneities. However, one is not free to arbitrarily specify the large scale effects 
of backreaction: They must be shown to arise from averaging an actual inhomogeneous 
spacetime. 

Our approach — which we summarize in detail in section Sec. |H] — bears considerable re- 
semblance to that of Zalaletdinov [6j in that we define a backreaction tensor fi a bcdef, which 
is closely related to Zalaletdinov's Z a b c ^ e f. However, we make full use of Einstein's equation, 
and we have better mathematical control over the perturbations through our introduction of 
one-parameter families, which allow us to rigorously take a short wavelength limit. "Averag- 
ing" of nonlinear terms in Einstein's equation may then be rigorously defined via the use of 
weak limits. We were thus able to derive strong constraints on fi a bcdef, and we proved the fol- 



lowing theorem about dynamical behavior in the limit of short wavelength perturbations: If 
the matter stress-energy tensor satisfies the weak energy condition (i.e., has positive energy 
density, as measured by any timelike observer), then the leading order effect of small scale 
inhomogeneities on the background metric is to produce an effective stress-energy tensor 
that is traceless and itself satisfies the weak energy condition. In particular, backreaction 
produced by small scale inhomogeneities cannot mimic a cosmological constant. 

However, our work merely assumed one-parameter families with the properties stated 
in Sec. [TT] below; we did not prove existence of a family with nontrivial backreaction, thus 
leaving open the possibility that our assumptions are self-consistent only in the case of no 
backreaction. Objections also have been raised that the interpretation of the parameter 
appearing in our one-parameter families is unclear [l0|, and that our averaging scheme 
is "ultralocal" , and thus cannot capture effects due to matter inhomogeneities over finite 
regions 



11] . Therefore, it seems clear that it would be useful for us to provide a concrete 
example of a one-parameter family of spacetimes with nontrivial backreaction that satisfies 
our assumptions. Such an example would serve the multiple purposes of proving the self- 
consistency of our assumptions, illustrating the meaning of the parameter appearing in our 
one-parameter families, and showing explicitly how our theorems have direct implications 
for perturbations of finite amplitude and wavelength. 

The purpose of this paper is to provide two examples of one-parameter families of space- 



times with nontrivial backreaction. The first, analyzed in Sec. IHIt is adopted from 12] and 
is a family of vacuum Gowdy spacetimes. This family satisfies all of the assumptions of our 
work and has nontrivial backreaction in that the limiting metric does not satisfy the vacuum 
Einstein equation. Its effective stress-energy tensor can be explicitly seen to be traceless and 
satisfy the weak energy condition, as guaranteed by the theorems in [3] . 

The second family, which we produce in Sec. [TV] is obtained by "Synge's method," i.e., 
by simply writing down a family of metrics and declaring the matter stress-energy of each 
spacetime to be equal to the Einstein tensor of the metric, thereby trivially "solving" Ein- 
stein's equation [13]. This family is non- vacuum and approaches a FLRW spacetime as the 
parameter approaches its limiting value. It satisfies all of our assumptions except that the 
matter stress-energy does not satisfy the weak energy condition. We obtain a nontrivial 
effective stress-energy tensor, which does not satisfy the properties of our theorems. This 
example shows that if matter violates the weak energy condition (or if Einstein's equation is 



not sufficiently utilized), then it is easy to produce spacetimes with significant backreaction 
that is unconstrained by our theorems. This illustrates the necessity of applying the full 
content of Einstein's equation with physically reasonable energy conditions imposed upon 
matter in any analysis of backreaction. 

In this work, we follow all notation and sign conventions of [14| . 

II. FRAMEWORK 

In this section we review our framework and we summarize the main results of 3j. Our 
ramework is a generalization to the non- vacuum case of a framework proposed by Burnett 



15] , which itself is a rigorous version of Isaacson's approach [161 ]. 

As indicated in the Introduction, to analyze a spacetime with backreaction g a b{^o), we 
attach to it a one-paramater family of spacetimes 5w,(A) for < A < A . We require g a b{X) 
to be a solution to Einstein's equation for all A > 0, and, as A — > 0, we require uniform 
convergence of g a b(X) to a "background metric" g ab = g a &(0). One-parameter families of 



metrics are also emp 



oyed in this manner to properly treat ordinary perturbation theory 



(see, e.g., Sec. 7.5 of 14(). However, whereas in ordinary perturbation theory the metric is 
required to be jointly smooth in A and the coordinates x, here we require only that spacetime 
derivatives of (g a b(X)—g ab ) are bounded as A — > 0. This elevates the importance of derivatives 
of perturbations of the background metric, so that a priori, small scale inhomogeneities 
can play a dynamical role in the evolution of the background metric. Indeed, in contrast 
to ordinary perturbation theory, it is not true, in general, that the Einstein tensor has a 
uniform limit as A — > 0, and it does not follow that g ab is a solution to the Einstein equation. 
However, if we add a suitable assumption that "spacetime averages exist" , then we can derive 
an equation for g ab , which contains the contributions arising from the backreaction produced 
by small scale inhomogeneities. 

Our precise assumptions are as follows. Fix a spacetime manifold M with derivative 
operator V a , and let g a b{X) be a one-parameter family of metrics on M, defined for A > 0. Let 
tab be an arbitrary Riemannian metric on M, and define |t a i...a„ | 2 = e ai&1 . . . e a " bn t ai ... an tf )1 ...b, i . 
Suppose now that the following assumptions hold: 

(i) For all A > 0, we have 

G ab (g(X)) + Ag ab (X) = 8vrT a6 (A), (2.1) 



where T ab (X) satisfies the weak energy condition, i.e., for all A > we have T a b(X)t a (X)t b (X) > 
for all vectors t a (X) that are timelike with respect to g a b{X). 

(ii) There exists a smooth positive function C\(x) on M such that 

IMA,aO|<ACi(aO, ( 2 -2) 

where h ab (X, x) = g ab (X, x) - g ab (0, x). 

(iii) There exists a smooth positive function (^(a?) on M such that 

\V c h ab {X,x)\<C 2 {x). (2.3) 

(iv) There exists a smooth tensor field Habcdef on M such that 

w-lim [V a h cd (X)V b h ef (X)} = /J abcd ef, (2.4) 

A — ^U 

where "w-lim" denotes the weak limit. 

The notion of "weak limit", which appears in assumption fliv| . corresponds roughly to 
taking a local spacetime average, and then taking the limit as A — > 0. More precisely, 
Ah...o„(A) converges weakly to A ai ...a n as A — > if and only if, for all smooth tensor fields 
yai...a n Q f com p ac t support, 

KmJ r- a "A ai ... an (X) = J r- a "A a °l an . (2.5) 

Assumptions (p|)-( pv|) allow us to derive an equation satisfieqj by g ab . Indeed, in |3|] we 
showed that g a b(0) satisfies the equation 

GM^ + ^i^SnT^ + Sntit (2.6) 

Here, T ab = w-lim a-s-o ^(A) — which must exist as a result of the assumptions — and t ab is a 
particular linear combination of contractions of the tensor [labcdef- The tensor rv is called the 
"effective gravitational stress-energy tensor," and it describes the dominanto contribution 
to backreaction due to small scale inhomogeneities. 
In [3J, we proved two theorems constraining t ab . 



1 Note that assumption Q only requires that the Einstein equation hold for A > 0, but not for A = 0. 

2 Higher order contributions to back reaction were also derived in [3( . In [17[ , we showed that the dominant 

effect of the higher order contributions on the dynamics of our universe is to modify the Friedmann 
equations by inclusion of effective stress-energy contributions associated with Newtonian gravitational 
potential energy and stresses and with the kinetic motion of matter. 
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Theorem 1 Given a one-parameter family g a b(X) satisfying assumptions (R|)-(lw|) above, the 
effective stress-energy tensor t ab appearing in equation (12. 6 p for the background metric g ab 
is traceless, 

t (0)a a = 0. (2.7) 

Theorem 2 Given a one-parameter family g a b(A) satisfying assumptions B\)- Kiv}) above, the 
effective stress- energy tensor t ab appearing in equation (12. 6p for the background metric g ab 
satisfies the weak energy condition, i.e., 

t^t a t h > (2.8) 

for all t a that are timelike with respect to g ab . 



In essence, these theorems show that only those small scale metric inhomogeneities corre- 
sponding to gravitational radiation can have a significant backreaction effecto In particular, 
in the case where we have FLRW symmetry, the effective stress-energy tensor must be of 
the form of a P = |p fluid, and therefore cannot mimic dark energy. 

Finally, we note that in the statement of assumptions fp])-f lrv|) above, A is a "continuous 
parameter," i.e., it takes all real values in an interval [0, Ao] for some Ao > 0. However, there 
would be no essential change if we took A to be a "discrete parameter" in our assumptions — 
such as A = {1/N} for all positive integers iV — provided only that A takes positive values 
that limit to 0. It will be convenient to use such a discrete parameter in the example of the 
next section. 

III. BACKREACTION IN VACUUM GOWDY SPACETIMES 

Gowdy spacetimes are exact plane-symmetric vacuum solutions which describe closed 
cosmologies containing gravitational waves 19j. We restrict to the case where spatial slices 



have topology T 3 . The general T 3 Gowdy metric may be written in the form 20] 



^4owd y = e (T ~ a)/2 {-e- 2T dr 2 + d0 2 ) +e~ T [e p da 2 + 2e p Qdad5 + (e p Q 2 + e~ p ) dd 2 ] , (3.1) 



3 Although matter inhomogeneities can produce only a negligibly small effect on the dynamical evolution 
of the background FLRW metric, they can produce large observable effects on the apparent luminosity of 
objects in the universe via gravitational lensing; see, e.g., 18 1. 



where the functions a, P and Q depend only on r and 9 and the spatial coordinates have the 
range < 9, a, 5 < 2ir, with periodic boundary conditions. We will restrict consideration to 
the case Q = 0, known as the polarized Gowdy spacetimes. 

When Q = 0, the vacuum Einstein equations reduce to the evolution equation, 

= P - e~ 2T P". (3.2) 

and the constraints, 

a = P 2 + e~ 2T (P') 2 , (3.3) 

a' = 2P'P. (3.4) 

Here, the prime and dot denote differentiation with respect to 9 and r, respectively. The 
general solution to ( 13. 2 p is 

oo 

P = A) + B r + J2 [ A nJo {ne- T ) + B n Y (ne~ T )] sin (n9 + <f> n ) , (3.5) 

n=l 

where A n , B n , and <fr n are free parameters. The discrete index n arises from the compactness 
of the 9 direction. 



Following 12J (see also 2l|). we now construct a one-parameter family satisfying assump- 
tions (TIT) (llvj) . As discussed in the last paragraph of the previous section, our family will be 
parametrized by a discrete parameter N, and the limit N — > oo will correspond to A — > 0. 
We choose 

P N = -=Jq (Ne~ T ) sin(A^e), (3.6) 

y N 



where A is an arbitrary constant. The constraint equations may be solved by direct inte- 
gration to yield 

A2 -T 

a N = -—J 1 (Ne- T )J (Ne- T )cos( y 2N9) (3.7) 



2?\r^-2T 



A 2 Ne 



{ [J {Ne-)} 2 + 2 [J, {Ne~ T )} 2 - J (Ne^) J 2 (Ne~ T ) } . (3.8) 



4 
From the asymptotic form of the Bessel function, we find that for large values of N, 



p n -> ^ V V cos { Ne ~ T ~ l) sin ( m ) ■ ( 3 - 9 ) 
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From this and the similar asymptotic form of ajy, it may be verified that assumptions (U|)- (!iv|) 
of the previous section hold for the family {P/v, a^}. In particular, as N — > oo, we have 

lim P N = 0, (3.10) 

N— >oo 

lim on = (3-H) 

N— >oo 7T 

where the convergence is uniform on compact sets. Thus, the "background metric" is 

ds 2 {0) = e (-+^- T /-)/2 (_ e -2r rfr 2 + d9 2j + e -r (^2 + ^2) (3^3) 



Despite the fact that for each N < oo, the metric is a vacuum solution, the limiting 
metric (I3.12p is not a vacuum solution. The effective stress-energy, r ab , could be computed 
by computing [i a bcdef via (12.41) and using the formula for t a1 j in terms of [i a bcdef given in 3j . 
However, it is far easier to obtain t ab by simply computing the left side of (12. 6p . It is easily 
seen that the nonvanishing components of the Einstein tensor of the metric (13. 1 2[) are given 
by 

A 2 e- T 



G TT 0? (o) ) 



An 
A 2 e T 



Gee (S (0) ) = -^- (3.13) 

Note that the effective stress-energy r ab = G a b/8n associated with ( 13. 13[) is traceless and 
satisfies the weak energy condition, in accord with our general theorems. 

Thus, we have provided an explicit example of a one-parameter family wherein the small 
scale inhomogeneities produce a nontrivial backreaction. One can see explicitly in this 
example how, in our framework, the vacuum, inhomogeneous spacetime metric g a b{N) can 
be well approximated for large but finite iV by the non-vacuum, homogeneous spacetime 
metric (13~T2D . 

IV. BACKREACTION WITHOUT ENERGY CONDITIONS 

In this section, we provide an example of a one-parameter family that satisfies assump- 
tions (Jn])-([iy]) of Sec. [IT]but does not satisfy the requirement of assumption @ that the matter 
stress-energy tensor satisfy the weak energy condition. This illustrates the importance of 
the weak energy condition on matter for the validity of theorems [1] and [2] of Sec. [Til 

11 



For our one-parameter family, g a &(A), we choose metrics conformally related to a spatially 



flat FLRW metric g ab , i.e. 



(0) 



g ab (X) = &(\)g™ = n\\)a\r) Vl 



ab • 



where a(r) is chosen arbitrarily. We choose the conformal factor to be 

logfi(A) = XA 



:mi x) +sin (D +sin (xAi 



(4.1) 



(4.2) 



with A constant. Of course, 5W>(A) does not satisfy Einstein's equation with any known form 
of matter. Nevertheless, we may simply declare the existence of a new form of matter with 
stress-energy tensor given by 

T ab (X) = ^-G ab (g(\)) (4.3) 

for all A > 0. Then g a fc(A) is a solution of Einstein's equation (with vanishing cosmological 
constant, A = 0) for all A > 0. This procedure for "solving" Einstein's equation is usually 
referred to as "Synge's method" [l3j. It can then be easily verified that our one-parameter 
family ( 14. ip and (14. 2p satisfies all of the assumptions of Sec. [Til except that the stress-energy 
tensor of our new form of matter does not satisfy the weak energy condition for A > 0. 

By (12. 6p . the effective stress-energy tensor produced by the small scale inhomogeneities 
is 

1 



t 



(0) 
ab 



8n 



Gab(g 



(oh 



T, 



(0) 



ali 



-L w -\im[G ab (g^) - G ab (g(\))} 

Q7T A->0 



— w-lim 

8tt a->o 



2 V a V 6 log fi - 2^y°) cd V c V d log n 



(0) 

,b ' 

(0) (0)cd 



-2 (V a log O) (V 6 log D.) - gS/gW* ( V c log £1) (V, log Q) 



(4.4) 



Here, on the last line, we used the expression relating the two Ricci tensors (see appendix 

d of y ), 

Rab(g(\)) = Rab(g {0) ) - 2V a V b logfi - g^ g^ cd V C V d \ogQ 

+ 2 (V„ log Q) (V 6 log n) - 2g®g® cd ( V c log Q) ( V d log Q) , (4.5) 



(o) 



where V a is the derivative operator associated with g y ab ' . A straightforward calculat 



tiortl 



The only nontrivial weak limits involved are of the form w-lim,\->o cos 2 (x/A) = 1/2. 
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yields 

4? = ^ (4-6) 

4 0) = (4.7) 

Thus, our effective stress-energy tensor produced by small-scale inhomogeneities corre- 
sponds to a P = — |p fiuid. This faild_| to be traceless and fails to satisfy the weak energy 
condition, in violation of both theorems Q] and [2] of Sec. [Til This should come as no surprise, 
since we have failed to satisfy assumption ([TJ). This example thus shows that assumption (jTJ 
is essentially needed for the validity of our main results. 



We note that several analyses [7|, |22j within the general framework of [6j have found an 
effective P = — |p fluiqj. As in our example above, this result is inconsistent with the 
tracelessness of t ab (theorem CD of Sec. [UJ. We believe that the violation of the tracelessness 
of t ab in [7J, |22] is of a similar origin: As discussed in the Introduction, the full content of 
Einstein's equation has not been utilized, and the requirement that matter satisfy the weak 
energy condition has not been imposed. 
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